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The direction of outflows from filaments: 
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ABSTRACT 



Context. It is generally presumed that the outflows from a YSO are directed close to its rotation axis (i.e. along its angular momentum 
vector and orthogonal to any attendant accretion disc). Many YSOs are formed from dense prestellar cores embedded in filaments, 
and therefore the relative orientations of outflows and filaments may place a useful constraint on the dynamics of core formation. 
Aims. We explore this possibility, from the viewpoint of what it may tell us about the angular momentum delivered to a core 
forming in a filament. We stress that we are not here addressing the issue of the relationship of filaments and outflows to the 
prevailing magnetic field direction, although this is evidently also an interesting issue. 

Methods. We use data from the literature and the SCUBA archive to estimate the projected angles between 45 observed outflows and 
the filaments which appear to contain their driving sources. The distribution of these angles is then compared with model predictions, 
so as to obtain a statistical constraint on the distribution of intrinsic angles between outflows and filaments. 

Results. Using the limited data available, and neglecting any possible selection effects or correlations between nearby outflows, we 
infer that the observed outflows have a tendency to be orthogonal to the filaments that contain their driving sources. Specfically, in the 
cases where the directions of the filaments and outflows are well defined, we infer statistically that 72 % of outflows are within 45° of 
being orthogonal to the filament, and only 28 % are within 45° of being parallel to the filament. 

Conclusions. This suggests that the prestellar cores which spawned the YSOs driving the observed outflows had angular momenta 
which were approximately orthogonal to the filaments out of which the cores formed. We briefly discuss the implications of this for 
two proposed core formation mechanisms. 

Key words, star formation - outflows - molecular clouds - interstellar gas dynamics 



1. Introduction 



The substructure within a star-forming molecular cloud of- 
ten appears to be filamentary, and the dense prestellar cores 
out of which protostars condense are often embedded within 
filaments (e.g. Schneider & Elmegreen 1979; Hatchell et al. 
2005; Johnstone & Bally 2006; Kirk, Johnstone & Tafalla 
2007; Nutter & Ward-Thompson 2007; Kirk, Ward-Thompson 
& Nutter 2007; Muench et al. 2007; Goldsmith et al. 2008; 
Narayanan et al. 2008). The implication is that the growth of 
prestellar cores is fed mainly by material flowing in along fila- 
ments, and - although Hatchell et al. (2005) caution against 
it - this interpretation is supported by many numerical sim- 
ulations, (e.g. Passot, Vazquez-Semadeni & Pouquet 1995; 
Nagai, Inutsuka & Miyama 1998; Padoan & Nordlund 1999; 
Balsara, Ward-Thompson & Crutcher 2001; Klein, Fisher & 
McKee 2001; Gammie et al. 2003; Li, P. et al. 2004; Li, Z. & 
Nakamura 2004; Nakamura & Li 2005; Oishi & Mac Low 
2006; Ciolek & Basu 2006; Kudoh et al. 2007; Kudoh & Basu 
2008; Hennebelle et al. 2008; Offner, Klein & McKee 2008). 
The highly anisotropic inflow from a filament onto a core may 
then have consequences for the dynamics of core collapse and 
fragmentation. In particular, the net angular momentum of the 
inflowing material will strongly influence the orientation of bi- 
nary orbits and circumstellar accretion discs in the small-Af sub- 
cluster of protostars forming in a prestellar core. 



During the embedded Class and Class I phases of proto- 
stellar evolution (e.g. Di Francesco et al. 2007; Ward-Thompson 
et al. 2007), the orientations of circumstellar accretion discs can 
be inferred from the directions of the outflows which they drive. 
An outflow is presumed to be driven by torsional MHD waves 
propagating along magnetic field lines anchored in a circumstel- 
lar disc (e.g. Pudritz & Norman 1986). Therefore the outflow 
should be approximately parallel to the rotation axis of the cir- 
cumstellar disc, i.e. along the angular momentum vector of the 
material forming the parental core. 

It follows that the orientation of an outflow, relative to the 
filament in which its driving source has been born, may con- 
tain important information on the dynamics of core formation. 
Specifically it constrains the relationship between the inflow 
forming the core (if we assume that this inflow is concentrated 
along the filament) and the angular momentum which this ma- 
terial brings with it (if we assume that this angular momentum 
is oriented along the direction of the outflow). This information 
can be used to discriminate between different mechanisms for 
core formation (e.g. Whitworth et al. 1995; Banerjee, Pudritz & 
Anderson 2006; Banerjee & Pudritz 2007). 

The present paper is organised as follows. In Section 2, we 
present the distribution of projected angles, y, between observed 
filaments and the outflows from YSOs that appear to be embed- 
ded in the filaments. We show that the distribution of projected 
angles is compatible with the distribution expected if outflows 
are approximately orthogonal to filaments; relative to the orthog- 
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Fig. 1. The open circles give the cummulative distribution of the 
observed (i.e. projected) angles y between filaments and out- 
flows. The full line gives the mean cummulative distribution for 
the best-fitting model. This model has cr^ — 0.8, where iff is the 
intrinsic angle between the outflow and the normal to the fila- 
ment. The dashed lines give the +lcr dispersion about this line. 



onal direction, the intrinsic directions of outflows have a stan- 
dard deviation <x^ ~ 45°, which implies that 72 % of outflows 
are within 45° of being orthogonal to the filament which contains 
their driving source. In Section 3, we discuss two mechanisms 
for core formation, (i) gravitational fragmentation of a shock- 
compressed layer (e.g. Whitworth et al. 1995), and (ii) gravo- 
turbulent fragmentation (e.g. Klein, Fisher & McKee 2001; 
Banerjee, Pudritz & Anderson 2006; Banerjee & Pudritz 2007; 
Offner, Klein & McKee 2008). We infer that the first mech- 
anism is more likely to produce the observed distribution, but 
this is not a very strong inference, and more observational 
data are required to make it more robust. We also discuss the 
statistical significance of the result, and the possible influence 
of selection effects. The data used are presented in Appendix 
A. The distribution of projected angles expected is derived in 
Appendix B. 

2. The distribution of observed and intrinsic angles 
between filaments and outflows 

On Figures \T\ and [2] open circles show the cummulative dis- 
tribution of y, the observed (i.e. projected) angle between out- 
flows from YSOs that appear to be embedded in filaments and 
the filaments themselves. The sources of the raw data, and the 
techniques used to derive the angles, are described in detail in 
AppendixlAl 

To interpret this data, we assume (i) that the intrinsic angle, 
iff, between the outflow and the normal to the filament has a trun- 
cated Gaussian distribution with standard deviation cr^ (trun- 
cated in the sense < < n/2); and (ii) that the observed 
systems are randomly oriented relative to the observer's line of 
sight. Using this model we have performed, for different values 
of cr^, a Monte Carlo experiment to generate 20,000 independent 
random samples, with each sample comprising N — 45 artificial 
observed systems. Using these random samples, we determine 
the mean projected angle of the « th member of a sequence, fiy(ri), 
and its standard deviation, <r r («). 

On Fig. Q] the full line shows the cummulative distribu- 
tion for /i y (n), and the dotted lines show the cummulative dis- 
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Fig. 2. The open circles give the cummulative distribution of the 
observed (i.e. projected) angles y between filaments and out- 
flows. The full line gives the mean cummulative distribution for 
the best-fitting model (cr^ = 0.8). The dashed line gives the dis- 
tribution when all outflows are exactly orthogonal to the filament 
(i.e. o~ tp = 0). The dot-dash line gives the distribution when all 
values of iff are equally likely (i.e. cr^ = oo). The dotted line 
gives the distribution when the angle between the outflow and 
the filament (i.e. if/' = n/2 - iff) has a Gaussian distribution with 
cr r = 0.8. 
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0.5 1 1.5 
ORTHOGONAL TO FILAMENT 

Fig. 3. Polar diagrams illustrating the distribution of outflow di- 
rections (in one quadrant) for the three non-singular cases il- 
lustrated on Fig. [2] The line types are the same, i.e. full for 
cr^ = 0.8, dot-dash for cr^ = oo, and dotted for oy = 0.8 (where 
if/ = tt/2 - iff). For o> = 0.8, 72% of outflows are within 45° of 
the orthogonal to the filament. 



tributions for /i y («) + <x y (n), for the best-fit model, with cr^ = 
0.8 radians ^ 45°. We have identified the best fit by minimizing 
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where y(ri) is the n th observed angle; for the best fit, Q ss 0.28 . 
cr^ = 0.8 implies that there is a 72% chance that the outflow 
is within 45° of being orthogonal to the filament, and hence a 
28% chance that the outflow is within 45° of being parallel to 
the filament. 

On Fig. [2] we show the model distributions (/i y (n) only) for 
cfy = (all outflows exactly orthogonal to the filament; dashed 
line), cfy = 0.8 (best fit; full line), and tfy = oo (no prefered 
orientation; dash-dot line). In addition we show the results when 
the angle between the outflow and the filament, i.e. iff' = n/2-iJ/, 
has a Gaussian distribution with cfy< = 0.8 (dotted line). The 
very large disparity between this last model distribution and the 
observed one reinforces our conclusion that outflows are very 
unlikely to be even approximately aligned with filaments. 

Fig. |3] is a polar diagram illustrating the distribution of out- 
flow directions for the three non-singular cases illustrated on Fig. 
[2] The same line types are used, i.e. full for cr^ - 0.8, dot-dash 
for cr^ = oo, and dotted for oy = 0.8 (where if/ = n/2 - ifi). 

3. Discussion 

3.1. Selection Effects 

Our conclusion that outflows tend to be approximately orthogo- 
nal to the filaments spawning their driving sources may be cor- 
rupted by by two factors. 

First, the observed systems (filament, YSO and outflow) are 
not just small in number, but are also located in a small number 
of local star formation regions. The orientation of systems in the 
same star formation region may well be correlated, for example, 
due to a preferential direction of external compression, or due 
to a prevailing, approximately uniform, large-scale background 
magnetic field. Our assumption that the observer's line of sight is 
randomly oriented relative to each system is then not valid. This 
problem can probably only be alleviated by extending the obser- 
vational database. Indeed, Strom et al. (1986) and Duchene & 
Menard (2004) find that in the Taurus-Auriga star formation 
region the outflows - but not necessarily all the protostellar 
discs - tend to be parallel to the large-scale magnetic field. 

Second, it is possible that outflows that propagate approx- 
imately orthogonal to a filament are more - or less - readily 
identified than those that propagate approximately parallel to a 
filament. It is hard to argue rigorously one way or the other. Does 
the higher density encountered by an outflow propagating ap- 
proximately parallel to a filament inhibit its advancement and 
thereby make it harder to identify, or does it increase its visibil- 
ity by supplying the dense gas required for a dissipative working 
surface (i.e. an HH Object)? We are inclined to believe the lat- 
ter, in which case our conclusion should be reinforced, i.e. more 
than 84 % of outflows should be within 45° of the orthogonal to 
the filament. 

3.2. Implications for core formation 

In order to speculate on how the apparent preferential orientation 
of outflows approximately orthogonal to filaments might be in- 
terpreted, we consider two possible core formation mechanisms. 
These two mechanisms are not chosen because they are the only 
ones that have been proposed, but because, as far as we know, 
they are the only ones which make specific predictions concern- 
ing the relative orientation of filaments and discs (hence out- 
flows). 

In the first mechanism (Whitworth et al. 1995), two interstel- 
lar gas flows collide and produce a shock compressed layer. In 



the centre of mass frame of the layer, the collision is unlikely 
to be head-on, so the collisions partners have net orbital angular 
momentum, H, which is orthogonal to the collision velocity, v. 
Hence the layer tumbles about h = H/|H|. If the post-shock gas 
is cooler than the pre-shock gas, the layer eventually becomes 
gravitationally unstable, and breaks up, first into a network of fil- 
amets, and then into prestellar cores along the filaments. Because 
the layer is tumbling about h, the filaments are also tumbling 
about h. Consequently the streamlines delivering material to the 
core from the two directions of the filament become increasingly 
offset relative to one another (see Figs. 8 and 9 in Whitworth et 
al. (1995)) and the prestellar core acquires angular momentum 
approximately parallel to h and orthogonal to the filament. In 
this case, we might expect the outflows from YSOs formed in 
the core to be approximately orthogonal to the filament (appar- 
ently as observed). 

In the second mechanism (e.g. Klein, Fisher & McKee 
2001; Banerjee, Pudritz & Anderson 2006; Banerjee & Pudritz 
2007; Offner, Klein & McKee 2008), filaments are created in 
the shear flow between two colliding turbulent streams. In the 
purely hydrodynamic cases presented by Klein et al. (2001) and 
Banerjee et al. (2006), and in the magneto-hydrodynamic case 
presented by Banerjee & Pudritz (2007), the filament is aligned 
along the orbital angular momentum of the collision, i.e. along 
h. Consequently the filament tends to rotate about its long axis 
(rather than tumbling). The material flowing into a core from 
the two directions of the filament now delivers angular mo- 
mentum which is aligned approximately parallel to the filament. 
Therefore in this case we might expect the outflows from YSOs 
formed in the core to be approximately parallel to the filament 
(in apparent conflict with the observations). 



4. Conclusions 

We have analysed the limited number of protostellar systems in 
local star formation regions where a YSO embedded in a dis- 
cernable filament drives an outflow. In these cases, the distri- 
bution of angles between filaments and outflows appears to im- 
ply that outflows are usually approximately orthogonal to the 
filaments which contain (and presumably have fed) the cores 
which have spawned their driving YSOs. Given the chaotic 
and impulsive nature of multiple star formation in a prestel- 
lar core, one should expect some variance in the relative ori- 
entation of filaments and outflows. Specifically, it appears that 
72 % of outflows are within 45° of orthogonal to the filament 
which contains their driving YSO. This inference is not very 
robust, statistically, but, if confirmed, it would favour a model 
in which filaments were formed by gravitational fragmentation 
of shock compressed layers, as against one in which filaments 
were formed from shear instabilities in colliding hydrodynamic 
or magneto-hydrodynamic flows. Further observations of out- 
flows from filaments are needed to validate this conclusion. 
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Appendix A: The sample of observed filaments and 
outflows 

The systems we analyse are all in nearby star formation regions. 

First, outflows are identified using the catalogue of Herbig 
Haro Objects (HHOs) by Reipurth (1999). 

Second, for some of these HHOs, the driving YSO has been 
identified: in the Orion Integral Filament, by Chini et al. (1997) 



and Yu et al. (1997); in p Ophiuchus, by Allen et al. (2002); in 
Serpens, by Davis et al. (1999); in Taurus, by Mundt & Fried 
(1983), Strom et al. (1986) and Strom & Strom (1994); and in 
NGC1333, by Strom, Vrba & Strom (1976) and Davis et al. 
(2008). Knowing the positions of the HHO and the driving YSO, 
the direction of the outflow can be estimated. 

Third we have used dust continuum maps to determine 
whether the driving YSO is embedded in a larger filamentary 
structure. For this purpose we have used SCUBA 850//m data 
from Nutter & Ward-Thompson (2007) for the Orion Integral 
Filament, SCUBA 850 /mi data from the SCUBA Archive (Di 
Francesco et al. 2008) for p Ophiuchus and Serpens; SCUBA 
450 jum data from the SCUBA Archive (Di Francesco et al. 
2008) for NGC 1333; and SPITZER 160 /mi data from Nutter 
et al. (2008) for Taurus. 

In the cases where filamentary structure can be identified, 
we have determined the direction of the filament, in as objective 
a manner as possible, using the SEXTRACTOR clump-finding 
algorithm in STARLINK. This algorithm requires the user to 
specify a detection threshold, in units of the background noise 
(for example T = 6<x) and an isophotal radius (e.g. % - 2), 
which determines the degree of smoothing over neighbouring 
pixels. Maxima in the map are fitted with elliptical isophotes, 
and the semi-major axis of the ellipse then corresponds to the 
direction of the filament. We vary T and H until this direction is 
approximately idependent of the exact choice of 7~ and "R, and 
the direction corresponds reasonably with what the human eye 
sees. 

If 77 OUT is the direction of the outflow, and 77 FIL is the direction 
of the filament (both measured clockwise from north), then the 
projected angle between outflow and filament is given by 

y = MIN { |?7 OUT - ?7 FIL | , n - \r] om - T] HL \ } (A.l) 

In Table lA. ll we list, in columns 1 to 3, the name, right ascen- 
sion and declination for the HHOs used in our analyis; these are 
all from Reipurth (1999). In columns 4 to 7, we give the name, 
right ascension and declination for the driving YSO, plus the ref- 
erence for the identification of the driving YSO. In columns 8 to 
10, we give the angles 77 OUT , ?7 fil and y. The systems (filament, 
YSO and HH object/outflow) are grouped in eight different star 
formation regions, and the source and wavelength of the contin- 
uum map used to define the filaments is given after the name of 
the star formation region. 

Fig. I A. 1 1 shows the continuum maps used to identify fila- 
ments; the long lines with arrows at either end indicate the de- 
rived directions of the underlying filaments. The circles indicate 
the positions of YSOs, and the crosses indicate the positions of 
HH Objects; the lines joining these are the inferred directions of 
the outflows. 



Appendix B: The distribution of projected angles 
between filaments and outflows 

Without loss of generality, we define Cartesian axes (x, y, z) with 
unit vectors (i, j, k), oriented so that (i) the filament is along the 
z-axis, i.e. it has a unit vector 

fl HL = k; (B.l) 

and (ii) the outflow is in the (y, z)-plane. If the outflow makes an 
angle i// with the y-axis, then the outflow has a unit vector 

"out = cos Wi + sin(^)fe. (B.2) 
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S. Anathpindika and A. P. Whitworth: The direction of outflows from filaments: constraints on core formation 




Fig.A.l. False-colour map of the 850 fim continuum emis- 
sion from part of the Orion Integral Filament (Nutter & Ward- 
Thompson, 2007), showing the positions of MMS2, MMS3, 
HH331, and the ellipse fitted to the local portion of the filament 
in order to estimate its direction. 

In order to compute the projected (i.e. observed) angle, y, 
between the filament and the outflow, we assume that the ob- 
server is located at large distance (compared with the size of the 
filament/outflow system) in the direction of the unit vector 



n OBS = sin(0) cos(0)i + sin(0) sin(0)j + cos(0)k. 



(B.3) 



We can construct a two-dimensional Cartesian co-ordinate 
system on the observer's sky using the unit vectors 



k x n 
|kxn Q 



sin(0)i + cos(</>)j , 



(B.4) 



■ cos(#) cos(0)i - cos(#) sin(0)j + sin(#)k ; 



the orientation of these axes is such that axb = n 0BS , so that 
(notionally) a is the abscissa, and b is the ordinate. 

On the observer's sky, the unit vector along the filament is 
then 



n 



(n FIL ■ a)a + (n FIL • b)b 



|(n FIL -a)a + (n FIL -b)b| 
and the unit vector along the outflow is 



ii' = 

OUT 



(n 0UT • a)a + (n 





b)b 



|(n 0UT • a)a + (n OUT • b)b| 

cos(ip) cos(0)a + (sin(i/f) sin(f9)-cos(</<) cos(6>) sin(0)}b 



f 1/2 (<P,0,4>) 



(B.5) 



where 

/OA, e, </>)-■ 



■ cos 2 (i/<) cos 2 (0) + cos 2 (i/0 cos 2 (f?) sin 2 (0) 
+sin 2 (i/r) sin 2 (f9)-2 cos(i//) sm(i//) cos(f9) sin(f9) sin( 



(B.6) 



Hence the projected angle, y, between the filament and the 
outflow, is given by 



cos(y) = n' • n' 

v 1 J FIL OUT 

- cos(i^) cos(0) sin(0) + sin((/») sin(f9) 

= / 1/2 OA,0,0) 



(B.7) 



We shall assume that iff is distributed according to the trun- 
cated Gaussian 



Pi, dip 



where 



C 



(2;r) I /V, 



C 



• - f 



exp 



1 



£_ 



exp^- 



dijj, 



m < (B.8) 



2<x 2 



dip . 



(B.9) 



If there are no selection effects, and none of the observed 
outflows are correlated in direction, then - again without loss 
of generality - putative observers can be distributed isotropically 
over one octant of the sky, with 



,d6d(f) = sm{6)d9 



2d</> 



< 6,i 



(B.10) 



For any given cfy (measuring the dispersion of the outflow 
direction about the orthogonal to the filament), the cumulative 
distribution of projected angles can be computed by means of a 
Monte Carlo integration. To do this we consider many indepen- 
dent sequences, each of N — 45 randomly distributed sample 
points in the integration space {ip, 6, <p). In this way we obtain 
both the expected cumulative distribution, P(y), and a measure 
of the dispersion about P(y) for a finite sample of N observa- 
tions. 

Specifically, we record the n th y-value from each sequence, 
and then compute their mean, fi y (n), and their standard devia- 
tion, cr y (n). The expected cumulative distribution, P(y), is then 

a plot of (n - IN against fJ. y (n). The lcr dispersion about 

this distribution is represented by the plot of (n - |) IN against 
fiy(n) + o-y(n). The distributions obtained in this way for various 
different cr^ are presented, and compared with the data, in Fig. 1. 



